Decay and storage of multiparticle entangled states of atoms in collective thermostat 
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We derive a master equation describing the collective decay of two-level atoms inside a single mode cavity 
in the dispersive limit. By considering atomic decay in the collective thermostat, we found a decoherence-free 
subspace of the multiparticle entangled states of the V^-like class. We present a scheme for writing and storing 
these states in collective thermostat. 

PACS numbers: 03.67.Pp, 03.65.Yz, 03.65.Ud 



O 
O 

> 
O 

o 



> 



o ■ 

^ : 

-i— > ■ 

G ■ 
cd 

3 : 



I. INTRODUCTION 

When information is encoded in a quantum state of the 
physical system, the robustness of the state is an important 
factor for successful communication. Due to decoherence, 
i.e., interaction with the environment, the state of the sys- 
tem can degrade and lose its quantum correlations. One of 
the possible solution of the decoherence problem is to use 
decoherence-free subspaces (DFSs); these include wave func- 
tions immune to decoherence . 

The first DFS has been introduced by Zanardi et al. |2[] 
for two-level atoms interacting with an electromagnetic field 
playing the role of environment. The wave functions belong- 
ing to the DFS are annihilated by the interaction Hamilto- 
nian and, therefore, are left invariant during evolution. Ex- 
amples of DFSs for various physical systems, in particular, 
for light, have been proposed by several authors (see, for ex- 
ample, GHH). 

Weinfurter et al. have demonstrated experi- 
mentally decoherence-free quantum communication based on 
the four-photon polarized states ((J. 

Simple observations show that quantum correlations be- 
tween particles can be produced and maintained in collec- 
tive processes. These are interesting for DFSs and a large 
number of physical systems with collective interactions can 
be found. For the Dicke model with a single resonant mode 
Bonifatio et al. have found a master equation describing 
a collective atomic decay when illuminated with a resonant 
mode [7]. Palma and Knight et al. |8[] have shown that two- 
atom decay can result in pure entangled states in a collective 
squeezed thermostat. However, entanglement of two atoms 
can be achieved in collective decay with vacuum thermostat, 
when atoms are placed inside a cavity, as it has been shown 
by Basharov [9]. There is a simple reason for collective decay 
in the cavity scheme. If atoms interact with a single cavity 
mode, then atomic relaxation arises because radiation leaves 
the cavity. This can be modeled as interaction between the 
mode and an external broadband field, which plays role of a 
thermostat. Therefore, the atoms, being coupled with the sin- 
gle mode, have a collective decay. 

The aim of this paper is to investigate the collective decay 
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of atoms in the entangled state. To achieve this we consider 
a simple model of two-level atoms inside a single mode cav- 
ity and a broadband field. For this model Klimov et al. ifloll 
derived master equations for atoms and a collective relaxation 
operator in the dispersive limit, assuming a vacuum bath. We 
discuss a more general squeezed thermostat, for which a mas- 
ter equation is derived using a formalism of unitary transfor- 
mations. There are some differences between relaxation oper- 
ators due to methods of derivation of master equations. From 
the physical point of view, in the dispersive limit there is only 
an exchange of phase between atoms and cavity mode. This 
is described by the effective Hamiltonian we found, which is 
diagonal over atomic and field variables. However, the behav- 
ior of atoms becomes more complicated due to an interaction 
between the cavity mode and broadband field, which plays the 
role of a thermostat. Then one finds a coupling of atoms with 
the thermostat, providing atomic collective relaxation. In this 
aspect our results differ from that of Ref. II 1 Oil . 

In contrast to ifloll we use the idea developed in Ref. ITlll ; 
we first find the effective Hamiltonian of the total system in- 
cluding the thermostat and then derive the master equation. 
Using the master equation we consider the dynamics of a class 
of multiparticle entangled states, which is a slightly general- 
ized W class introduced by Cirac et al. lfl2tl . Some of the op- 
tical and atomic implementations of the presented states have 
been demonstrated experimentally by Weinfurter et al. lfl3tl 
and Schmidt-Kaler 11411 . Some properties of these states, dif- 
ferent schemes to generate them, and several applications have 
been considered in Ref. lfl5ll . We also find that in the case 
when these entangled states are reduced to the Dicke states, 
they belong to DFS and are immune to collective decay. To 
explain this feature we use symmetry arguments. In fact, the 
total space of the Dicke states is represented by irreducible 
subspaces distinguished by their symmetry type. The collec- 
tive interaction we consider does not mix the wave functions 
from different subspaces due to symmetry conservation. Us- 
ing these properties, we present a model of a quantum mem- 
ory for writing, storing and reading information encoded in 
these entangled states. 

The paper is organized as follows. First we derive the mas- 
ter equation in the dispersive limit assuming a general model 
of thermostat. Then we introduce a set of multiparticle en- 
tangled states which can be reduced to the Dicke family and 
consider their decay in squeezed and vacuum thermostats. Fi- 
nally we present a scheme for reading and storing entangled 
states in the collective thermostat. 
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II. INITIAL EQUATIONS 

By considering the interaction between atoms and a field, 
one can obtain a master equation for one of the systems. This 
equation is also known as kinetic and often has Lindblad form. 
It describes irreversible processes including atomic relaxation, 
absorption or amplification of light, and others phenomena 
which can be reduced to the Lindblad equations. 



A. Hamiltonian 

We consider n two-level atoms inside a high-finesse optical 
cavity, a single cavity mode, and a broadband field outside the 
cavity. We assume the Hamiltonian of the system has the form 

H = H a + H c + H b + V 1 + V 2 , (1) 

where the Hamiltonians of free atoms, cavity mode, and 
broadband field are, respectively, H a = huioR^, H c = 
fwjc, H b = here R 3 = £ 3 (|1) 3 <1| - 

|0)j(0|), 0, 1 label the lower and upper levels of atom, and 
c, c\ b u , b^ are creation and annihilation operators for pho- 
tons of the cavity mode and broadband field, respectively. The 
interaction between atoms and cavity mode V\ has the form 

V 1 =g(jR- + cR+), (2) 

where collective atomic operators are given by R± = 

J2 3 R±\R { + = [R U) ] f = |l)j(0|. The term V 2 describes 
two processes: (1) an interaction between the broadband field 
and the cavity mode due to nonzero transmittance of the out- 
put mirror; (2) an interaction between atoms and the broad- 
band field due to non-ideal sidewalls of the cavity. It reads 

V 2 = X)M r "C f +J2 K "J R + ) ]+ h - c - ( 3 ) 

w 3 

From the physical point of view the broadband field plays the 
role of a thermostat and causes the relaxation of atoms and 
cavity mode. Relaxation terms can be achieved by switching 
on this field. It can be done in different ways using various 
approximations. 



B. Dispersive limit 

We assume that detuning A = \uj c — ujq\ is large and con- 
sider the dispersive limit, which can be justified when 1U7I1 

|A| » ngy/(ctc) + l. (4) 

To derive the master equation let us introduce a transforma- 
tion of Hamiltonian H given by a time independent unitary 
operator S 

H' = e~ iS He iS = -i[S; H] - (1/2) [S; [S; H}] + . . . . (5) 



Using perturbation theory over interactions V\ and V% one 
finds the operator S, from which an effective Hamiltonian de- 
scribing the interaction between atoms and cavity mode can 
be obtained. This Hamiltonian is diagonal over the field and 
atomic variables and has the form 

2 R-R+ + cJ2R 3 
H e = g ^ . (6) 

Under this approximation there is another effective Hamilto- 
nian H g which describes the interaction between atoms and 
broadband field 

H g = -^J2 r ^ R + b " + R - b ^- ™ 

In contrast to the second term in V 2 the obtained Hamilto- 
nian H g describes the collective interaction of atoms. In- 
deed, in the usual case of the dispersive limit there is no 
energy exchange between atoms and light, and this is in ac- 
cordance with the effective Hamiltonian H e , which is similar 
to Ref. [10fl. In the same time the effective Hamiltonian H g 
shows that atoms and the thermostat field exchange excita- 
tions. This is a particular feature of the dispersive limit due 
to the initial interaction ([3]). A close analogy is parametric 
down conversion in transparent nonlinear media, in which the 
virtual transitions result in an interaction between photons. 

Now we have a problem specified by H' = H a + H b + 
H c + H e + H g + V 2 , where broadband field can be considered 
as a thermostat in a given state. Assume the thermostat is 5- 
correlated and its state is squeezed with a center frequency f2: 

{blb u >) = N{u)5 u ^, 
(b^bl) = (N(u) + l)6 u , u ,, 
(b u b u >) = M(w)£ 2 n,a>-u>', 
(&t^,)=M>)<W+"'= (8) 

where the photon numbers N(u>) and M are related as 
\M(uj)\ < y/N(u)[N(u) + i\. A physical model of 
this thermostat can be represented by the light generated 
in parametric down conversion process. Its simple non- 
degenerate version is described by the next Hamiltonian H = 
^ u] {kLjbl 1+LU b n _ UJ + H.c), where 2il is the pump frequency 
and uj belongs to a frequency band h given by phase match- 
ing conditions. The photon numbers N and M have the form 
N{uj) — sinh 2 r w , M(u>) — exp(i argfc w ) coshr^ sinhr u , 
where r w ~ \k u \ is a squeezing parameter. For a squeezed 
vacuum r«l and N w 0, while M w exp(i arg fc w )r w . The 
generated light is broadband if h is much bigger than all repre- 
sentative frequencies of the problem, like the atomic and cav- 
ity mode decay rates. More precisely, assume that the width 
of the squeezed broadband field given by Eqs. © is much big- 
ger than the detuning A as in Then following the standard 
procedure of replacing a finite bandwidth system with white 
noise [ 16[] we can make all parameter of squeezed light to be 
independent from the frequency: N(u>) — N, M(u>) = M. 
We assume the squeezed thermostat is modeled by a paramet- 
ric down conversion source. Then its bandwidth is determined 
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by the phase matching conditions, which can be experimen- 
tally varied on a wide range. The next step is switching on 
the broadband field. This can be achieved by several methods 
based on projection operator techniques, stochastic differen- 
tial equations, and others. In any case we need a Markovian 
approximation to obtain a closed equation. In our case this 
means that the evolution of the broadband light is given by 
the free Hamiltonian Hi, only and the thermostat parameters 
TV and M are frequency independent. As a result we find a 
master equation for the density matrix p of atoms and cavity 
mode. The equation includes an effective Hamiltonian H e and 
relaxation terms. In the dispersive limit and in the interaction 
picture the master equation has the form 

p=-(i/h)[H e ;p}-Tp, (9) 

where the relaxation operator T includes three terms of the 
Lindblad form: T — £ • £j + C c + C a - The first term de- 
scribes the independent decay of atoms in the squeezed ther- 
mostat. When the atoms have the same coupling constant 
Ka,j = K u it reads 

C jP = { 1 j2){R ( ( ) R { l ) p - 2R i l ) pRf + pR^R^) 
+ ( 7T /2)(i? ' ) i?^ ) P - 2R<f ) pR U) + P rV>rV>) 
-2MK 2 R { i ] pR { l ] - 2M*K 2 R U) pR U \ 

where the decay rates of atomic levels are denoted by 7^ = 
|iq 2 (iV+l), 7T - \K\ 2 N,Bind\K\ 2 = »- 2 £jX w | 2 <J(w - 
uS), In free space one finds that \K\ 2 reduces to the well- 
known formula for the spontaneous decay rate Auo 2 d 2 / 3 He 3 . 
Equation (TlOb describes spontaneous decay of independent 
atoms in the squeezed thermostat, for which the transver- 
sal decay rate becomes slow because of squeezing: j± = 
(71 + 7t)/ 2 - Re{MK 2 }. The second term of T is the re- 
laxation of the cavity mode due to photons leaving the cavity, 
and has the form 

C c p — \T\ 2 [(N + l)(c^cp — 2cpc t + pc^c) 
+N(cc)p - 2c) pc + pec 1 ) 
+M(ccp — 2cpc + pec) 
M*(Jc*p - 2c t pc t + fxfc% (10) 

where |T| 2 = H~ 2 |r w | 2 5(w c — u). If R is a reflectance 
of the output cavity mirror, then |T| 2 — > c(l — R)/2L, where 
L is length of the cavity. The collective decay of atoms is 
represented by the operator £ a : 

£aP = \x?[{N + l)(R+R-p - 2R_pR + + pR+R_) 
+N(R-R+p - 2R+pR_ + pR_R+) 
+M(R+R+p - 2R+pR+ + pR+R+) 
+M*{R^R-p - 2R_pR_ + pR_R_)],(ll) 

where \x\ 2 — \gT/h 2 A\ 2 r, r = L/c. As a result, in the dis- 
persive limit there are three relaxation operators describing 
single-particle and collective decay. They have a straightfor- 
ward physical meaning and they differ from the relaxation op- 
erator in Ref. 1 1 Oil - which has cross terms including products 



of the collective atomic operators by operators of the cavity 
mode. 

In order to consider the collective decay of 
atoms let us introduce the interaction picture p 1 — 
cxp(—iftr 1 H e t)p exp(i?i _1 H ~ e t) and assume the following 
approximations: Let the first term in H e and single-particle 
relaxation be small, g 2 R-R + /HA, J2j £j ^ cc^R^/hA, 
C a , C c . TMsistraeif fl a /ft|A|, lu « (cx*)g 2 /\A\, \ X \ 2 n, 
and |r| 2 ( C t c )/n. Then we can neglect the difference between 
H e and g 2 2c [ cR^/hA so that the master equation for the 
atomic density matrix / = Ti c p' is 

/ = -£«/• (12) 



III. COLLECTIVE DECAY AND STORAGE OF 
ENTANGLED STATES 

When considering decay of atoms one finds that quantum 
correlations between particle can be supported in collective 
thermostats and the final or steady state depends from the ini- 
tial one. 



A. Entangled Dicke states 

Let's introduce the multiparticle entangled states, the slight 
modification of the W states discovered by Cirac H2I1 

t] n (l) = qi \10 . . .0) + fc|01 . . .0) + • • • + ?n |00 . . . 1),(13) 

where £ fe \qk\ 2 = 1- Some of these states are reduced to the 
Dicke states \ jma) lfl8ll . specified by three quantum numbers 
j, m, a, where |m| < j = 0, . . . , n/2 — 1, n/2, n is a number 
of particle and parameter a describes the degeneracy and takes 
rij = C?J 2+: ' — C^ / ' 2+J+1 values. The numbers j and m 
are eigenvalues of two commuting collective operators J3 and 

J 3 \jma) = m\jma), J 2 \jma) = j(j + l)\jma), (14) 

where Jf, obeys the commutation relations of the momentum 
operators [Jf,; J c ) = iebcdJd, b,c,d = 1, 2, 3. In the consid- 
ered case Ji = (1/2)(J2_ + R+),J 2 = (£/2)(ii+ - R-). 
When 

X> = 0, (15) 

k 

then there is a set of the zero sum amplitude states discovered 
by Pati IU9I1 , However, the wave functions rj n (l) under con- 
dition (fT3T > belong to the Dicke family with j=m=n/2-l ifsoTl . 
The states have the next representation 

n 

Vn(l) = V2^2q k \^-) lk ^\0) m , (16) 

k=2 

where |0)(im denotes a state of n — 2 particles (without first 
andfcth),*- = (1/V2)(|01) - |10)). Equation (O gives the 
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structure of entanglement of rj n (l); it tells that one of the par- 
ticles, say, 1 forms EPR pairs with all other particles 2, . . . , n 
and this feature is invariant under permutations of particles. 
Due to antisymmetric vectors the collective evolution of 
n particles in the state (TToT l involves only n — 2 particles. This 
has a simple reason. Two-particle collective operators R± and 
i?3 annihilate then for any operator of evolution U de- 
pending on R± , i? 3 we have 

U\ Vn )(r, n \ = 2^ |*- fc )(*rJ£/(lfc; ls)|0) (lfc)(ls) (0|, (17) 

ks 

where U(lk; ls)|0)n^(i s )(0| acts on all particles except for 
1, k and 1, s. These features allow us to get simple exact so- 
lutions for several problems of collective decay of rj n (l). 



B. Collective squeezed thermostat 

In squeezed thermostat there is an interesting feature. It 
can produce or store quantum correlations between particles 
for several initial states. Considering a two-atom collec- 
tive decay with initial density matrix /(0) = A|00)(00| + 
S|*+>(*+| + C|11)(00| +C*|00)(11| + D|11)(11|, where 
= (1/V2)(|01) + 110)), A + B + D = 1, one finds the 
next pure steady state H 

a = (VN + 1|00) + VJV|11))/(V2JV + 1). (18) 

This state is entangled. Note that this solution is correct for 
the initially symmetric state /(0). 

Consider the collective decay of the two entangled states 773 
and 774 described by (fT2l . Under conditions ( TT~5T > the wave 
functions read 

V3=q2*i 2 \0)3 + q3*i 3 \ h, 

m = 92*1210)23 + 93*1310)24 + ?4*ul°>23. (19) 

According to Eq. (T% the evolution of the density matrix 
1 773 ) (773 1 reduces to the dynamics of the single-particle state 
10) (0| for which there is a simple solution |0) (0| -> A|0) (0| + 
(1 - A)|l)(l|. The A = N/(2N + 1) is the occupation num- 
ber of the lower atomic level. One finds that 773 decays into 
a mixed state with complex structure. This fact can be ex- 
plained using a symmetry argument, which tells that under 
the single-particle decay the symmetry of the initial state is 
not conserved. In contrast dynamics of 773, the dynamics of 774 
has other features. To obtain the solution we use (TTBl and find 
that the final state is obtained by replacing 1 00) — > s 

V4 -> 92*r 2 l S )23 + 93*r 3 l S )24 + 94*r 4 l S )23- (20) 

From this equation it follows that the state is pure, has a more 
complicated entanglement structure, but as before one of the 
atoms forms EPR pairs with all another atoms. 



C. Vacuum thermostat 

Assuming a simple thermostat model for which M = N = 
0, the master equation ( TTZb reduces to 

/ = -K(R+R-f - R-JR+ + fR+R-), (21) 

where n — \x\ 2 - This equation describes a collective decay in 
the vacuum thermostat conserving quantum correlations. The 
simplest example is the two-particle antisymmetric function 
\f r_ belonging to DFS and immune for decay. The more in- 
teresting examples, introduced by Zanardi |2J], are DFSs of 
multiatom states, products of vp - . 

Suppose the atoms inside the cavity are prepared in the state 
77,1(1), then they evolve according to Eq. ED which can be 
solved exactly. It is easy to verify that the Lindblad opera- 
tor C f = R+R-f - R-FR+ + H.c. in Eq. (ED has the 
following properties: 

C \Vn)(Vn\ = Q|l;7i)<77„|-|Q| 2 |0)(0| + Q*|77 n )<l;n|, 
£ \r) n )(l;n\ = Q\l; n)(l; n\ - 2Qn|0)(0| + n\l; n)( Vn \, 
£o|l;n)<»k| = [£ |77„)(l;n|] t , (22) 

where Q — J^k |0) — |00 ... 0) is a ground state of atoms 
and 1 1; n) is a fully symmetric state, the normalized version 
of which, W n = (l/y/n)\l; n), is known as W state 

W n = (1/V»i)(|10 . . . 0) + |01 . . . 0) + . . . |00 . . . 1)). (23) 

It follows from Eqs. ( l22b that the Lindblad operator £ maps 
the set of states {|%)(t7„|, |0)(0|, |1; n){l; n\, |1;ti)(t7„|, 
|?7 n )(l;n|} into itself. This observation allows us to get an 
exact solutions for density matrix 

f(t)=A(t)\l;n)(r ]n \+A*(t)\r 1n )(l;n\ 
+B(f)|l;n)<l;n| + 5(t)|Q)<0| + D\r] n ){r) n \, (24) 

where the normalized condition reads A(t)Q* + A(t)*Q + 
B(t)n + S(t) + D(t) = 1 and coefficients obey equations 

A = -k(Au + DQ), 
B = -k{AQ* + A*Q) - 77iB, 
S = -2n(nS + n(l - D) + \Q\ 2 D), 

D = 0. 

Similarly to the squeezed thermostat there is a steady state 
solution, if t — ► 00 

/„ = D[-(Q/n)\l;n) + \r, n )][-(Q*/n)(l;n\ + ( Vn \] 

+ [(l-D) + \Q\ 2 D/n}\0)(0\(25) 

which depends on the initial state through parameter D. If 
D = 0, then f ss = |0)(0|. If D = 1 one finds evolution of 

Vn(l) 

\Vn)(Vn\ -» (Q/V^)(e- nKt - l)\W n )( Vn \ + h.C. 

+{\Q\ 2 /n)(l-e- nKt ) 2 \W n ){W n \ 
+ (|g| 2 /^(l - e- 2nKt ) \0)(0\ + \ Vn )( Vn \. (26) 
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It follows from (l26b that under condition (TT3T > T) n (l) is Dicke 
state and has immunity to the collective decay. However, this 
result can be obtained without any calculations because its an- 
nihilation by the Lindblad operator Co- In contrast to rj n (l), 
the fully symmetric W n state degrades: W n — ► |0). 

The robustness of r\ n entangled states can be clear from the 
symmetry argument. In the considered collective processes 
the particle permutation operator is an integral of motion, so 
that state symmetry is conserved. Therefore the antisymmet- 
ric wave function is robust to decay because the transi- 
tion — > |0) is forbidden, but the fully symmetric W state 
can transform into |0). In the case of r) n (l) the situation is 
more complicated, nevertheless symmetry plays a principal 
role here also. As it is known the space of Dicke states is rep- 
resented by irreducible subspaces distinguished by their sym- 
metry type over particle permutations. Under the condition 
(fT5T l the wave functions T] n (l) belong to [n, n — 1] irreducible 
representation of Dicke states in contrast to W and ground 
state |0), which belong to the [n, 0] one. Due to symmetry 
conservation the subspaces of different symmetry do not mix. 
This point is in accordance with the fact that in the dynamics 
of the wave functions the final states at t — > oo are not usual 
steady states but depend on their subspace and initial condi- 
tions. 

Robustness of r/„(l) is a natural basis for a quantum mem- 
ory. Memory includes writing, storing, and reading of in- 
formation encoded by a quantum state. By choosing rj n (l) 
states to encode information, a model of quantum memory 
can be designed. Writing and reading are achieved by swap- 
ping: a ® b — > b ® a. A particular case of swapping of two 
mode light in a Fock state into atomic ensemble has been con- 
sidered in Ref . 112 ill . Here we introduce a scheme for writing 
and storing multiparticle states. 

Assume that an interaction between atoms inside cavity and 
light, presented by its spacial modes with the wave vectors j, 
has the form 

V = ih^2 f( R + a j exp(ijrj) - h.c), (27) 

3 

where rj is a position of jth atom. This Hamiltonian describes 



an exchange of excitation between a single atom and a sin- 
gle mode. Assume that initially atoms and field are indepen- 
dent: |0) a ® \r] n )b, where |0) a is the ground state of atoms and 
| r] n )b is the light state given by Eq. d!31 l. where |0), |1) are the 
Fock states with and 1 photons, respectively. The multimode 
light can be prepared, i.e., by a set of beam-splitters, distribut- 
ing a single photon to different paths. For simplicity assume 
exp(ijrj) w 1, then the evolution is given by 

exp(-i^ 1 V r i)|0) Q ® \r) n ) b = cos(/t)|0) o <g> |?7 n ) 6 

+ sin(ft)\rj n ) a 0\O) b . (28) 

If sm(ft) = 1, the state of light \r) n )b is swapped into atoms. 
Under condition Sl5[ it can be stored in the collective ther- 
mostat. Due to the unitarity of transformation j28l , we can 
achieve a reading of the atomic state. 

IV. CONCLUSIONS 

Being collective properties of a physical system, quantum 
correlations between particles and entanglement can be pro- 
duced and stored in the collective processes. These processes 
can describe interaction between the physical system and its 
environment, which often plays a role of thermostat. In con- 
trast to the usual thermostat the collective thermostat supports 
quantum correlations and it is possible to find a DFS, which 
is a natural basis for quantum memory. For the considered ex- 
ample of collective decay atoms inside cavity, we found that a 
set of entangled states of the VF-like class is decoherence free 
and therefore is suitable to encode quantum information for 
storing it in collective thermostat. 
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